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Heated  From  Within 
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and  Yasuo  Shimazu 


Abstract 

In  the  present  paper,  the  Benard’s  cell  problem  in  a 
rotating  sphere  is  studied.  By  a preliminary  study  in  section  3, 
it  is  shown  that  the  rotation  stabilizes  the  convection  currents 
and  tends  to  make  the  scale  of  the  convective  cells  smaller.  In 
section  4 it  is  shown  that  each  poloidal  mode  of  the  motion  can, 
in  a non-rotating  sphere.-  be  excited  independently.  As  is  shown 
in  section  5,  this  is  not  the  case  In  a rotating  sphere.  Each 
motion  is  then  composed  of  infinite  numbers  of  poloidal  and 
toroidal  modes*  To  this  rather  complicated  case,  the  general 
results  in  section  3 can  still  be  applied.  In  connection  with 
the  present  problem,  a method  is  given  in  section  8 to  get  an 
equivalent  viscosity  coefficient  of  magnetic  fields. 


■&)  This  paper  (hitherto  unpublished)  was  completed  previous 
to  Dr.  Takeuchi's  joining  this  project.  It  is  here  reproduced  on 
account  of  its  importance  for  the  hydrom agne tic  theory  of  the 
earth’s  liquid  core.  (W.  M,  E.  ) 
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1.  In  a series  of  previous  papers  (H.  Takeuchi  and  Y.  Shimazu, 
1952,  1953).  it  was  shown  chat  a 3 elf -exciting  process  Is  possible 
by  which  the  earth's  main  magnetic  field  may  be  produced  and 
maintained.  The  self -exciting  process  is  considered  to  be  main- 
tained, by  the  Induction  currents  caused  by  the  motions  of  the 
fluid,  of  which  the  earth’s  core  is  composed.  In  order  to  make 

the  study  on  the  self -exciting  dynamo  complete,  however,  there 
remains  another  problem  to  be  solved.  That  is  the  problem  of  the 
fluid  motion  itself.  Is  it  possible  that  the  required  fluid 
motion  takes  place  in  the  earth's  core,  and  what  are  the  conditions 
for  this?  Is  the  fluid  motion  appropriate  to  maintained  the 
self -exciting  dynamo?  These  questions  will  be  considered  in 
the  present  paper. 

2,  After  considering  all  possibilities,  it  Is  now  believed 

that  the  fluid  motion  in  the  earth's  core  is  the  convective  one 

caused  by  non-homogenaous  heating  of  the  fluid  (E.  C.  Bullard, 

1949s  W.  M.  Elsasser,  1950).  In  fact,  our  previous  studies, 

referred  to  above,  are  based  on  this  convection-current  model. 

2C 

The  S2  -type  velocity  field  in  these  papers  is  nothing  but  the 
mathematical  expression  for  the  convection  current.  It  is  this 
SgJ~type  velocity  field  that  makes  our  self -excited  dynamo 
possible.  In  view  of  these  circumstances,  our  immediate  problem 
may  be  stated  as  follows: 

2 c 

(1)  Can'  a stationary  fluid  motion  of  -type  take  place 

in  the  earth's  core,  and  under  what  conditions? 

2 c 

(2)  It  must  be  shown  that  the  30  -type  velocity  field  is 

2c 

the  easiest  type  of  motion  to  excite..  Otherwise,  the  Sg  -type 


other,  more 


velocity  field  may  be  overcome  by  fluid  motions  of 
p 07.  e rf ul  t yoe  s . 

2 c 

The  stationary  velocity  field  of  S0  -type  reminds  us  of  the 
well-knov/n  Bsnard  cell.  It  may  be  that  the  fluid  motion  in  the 

. i 

earth* s core  is  of  a character  similar  to  the  Benard  cell  in 
a vessel  heated  from  below.  This  will  be  the  main  point  of 
view  of  the  present  study.  Many  studies  have  been  made  on  the 
convective  motion  of  a viscous  fluid  heated  from  below.  In  almost 
all  these  studies,  the  problem  is  treated  for  plane  boundaries, 
and  the  effects  of  rotation  upon  the  fluid  motion  are  net  taken 
into  account.  In  the  next  section,  taking  the  earth's  rotation 
into  account,  we  shall  treat  the  problem  for  plane  boundaries. 

3„  We  shall  consider  the  fluid  of  density  viscosity  ji, 

kinematic al  viscosity  v = thermal  diffusivity  -k  and  coeffi- 
cient of  thermal  expansion  a.  The  fluid  is  contained  in  two 
walls  placed  horizontally  at  a vertical  distance  of  . h • 
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The  fluid  is  heated  from  below.  The  temperature  gradient  in  a 
state  of  no  convention  is  p.  Referring  to  the  rotating  rect- 
angular coordinates  as  shown  in  Fig.  1,  we  have  the  following 
equations 

p §£  = - H+ 

Pm; = “ H + ^2v  " 2^u» 

pi* 


Qz. 


TVP  2 

^ T. 


(3.1) 


Customary  notations  are  used.  Attaching  the  suffixes  o to  the 
quantities  in  the  state  of  no  convection*  wo  have 


“o  = ’o  “ "o  = 0 * Po  - P„,z>  • To  “To<‘>* 


dT 


o _ 


dz 


= - pz  , /oo  = po>l  (1  - oT0), 


aP 


- f>  <S  ~ o, 
dz  1 o ° * 


(3.2) 


Pq  o ^eing  a certain  mean  density  corresponding  to  d certain  mean 
temperature  (or  height.,  z).  We  shall  put 


p - ?a  + p , T = v+  T * p =Pa  + /° 


(3.3) 


for  p,  T and  in  (3.1),  and  neglect  the  second  order  terms  in  u, 
v,  ws  p , T and  fo  . Furthermore,  we  shall  assume  that  the  density 
variation  of  the  fluid  is  caused  by  the  thermal  expansion  only  and 
that  the  fluid  behaves  otherwise  as  if  it  were  an  incompressible 
one.  Thus  we  get 

p v0i0d-«*i  . p » i3-4> 


Inserting  (3.2)  - (3.4)  into  (3.1)  and  denoting  p , T and^Q  Q 
as  p,  T and  p anew,  we  have  for  the  state  of  stationary  convection 

v^72u  - — ^ + 2 oj  v = 0, 

p ox  * 

vV2v  - ^ - 2<ou  = 0, 

■r—7  2 1 35 

vV  w - P + <v>  gT  = 0, 

k72T  + pw  = 0.  (3.6) 


While  the  above  equations  were  obtained  for  an  incompressible  fluid, 
H.  Jeffreys  (1930)  has  shown  that  the  same  equations  can  also  be 
applied  to  compressible  fluids  provided  the  density  does  not  vary 
greatly  within  the  system  and  provided  also  we  interpret  (3  as  the 
temperature  gradient  in  excess  of  the  adiabatic  one.  Putting 

(Us^^RT)  = e^mu  + n5r\u(z),  v(z),  w(z),  p(z),  T(z)^,  (3.7) 


we  get  five  ordinary  differential  equations  for  five  unknown 
functions  u(z),  v( z ) , - - - -.  Putting 

m2  + n2  = (irb)2  , z = ,b-  y (3.8) 


and  eliminating  u(z),  v(  z ) , w(z)  and  p(z)  from  (3.5)  - (3.7),  we 
have 


(«L  - ^2b2)3  + (2S&!)2  sL  + itm2i»Vne« 
llf  ’ v if 


T = 0.  (3.9) 


In  order  to  get  a quantitative  result,  we  shall  assume 

d2w 

T = w ~ = o 


(3.10) 


at  ^ = 0 and  1.  By  (3.10)  the  temperature  at  the  two  free  boui',dar ies. 


z - Cfi  and  z - h is  kept  constant.  These  boundary  conditions  are 
satisfied  by  putting 

T cC  sin  stc (s  = 1,  2,  — -).  (3.11)' 

Inserting  (3.11)  Into  (3.9),  ws  get  the  following  condition  for 
stationary  convection 


(b* 


2,3 

s ) 


+ ( 


2(gh^  ,2 

M / 


n 

3* 

7 


(3.12) 


where 


31 


kv 


(3.15) 


In  the  case  where  the  value  of  X in  (3.13)  is  larger  (or  smaller) 
than  that  given  by  the  right-hand  side  of  (3.12),  the  fluid  motion 
is  considered  to  be  unstable  (or  stable)  and  growing  (or  damped). 

— o • 

Keeping  ) as  a parametric  constant  and  making  use  of  (3.12), 

we  can  determine  the  critical  value  X as  a function  of  b and  s. 

The  value  of  X thus  considered  may  become  a minimum  for  certain  val- 
ues of  b and  s.  As  the  value  of  X in  (3.12)  i3  -a  monotonically 
increasing  function  of  3,  we  have  s = 1 as  the  value  of  s for 
which  X in  (3.12)  Is  a minimum.  Thus,  putting  a = 1 ans  ^ = 0 
in  (3.12),  we  can  determine  the  value  of  b for  which  X has  its 


minimum.  The  value  obtained  Is  an  increasing  function  of 
— 2 

( JL“~)20  In  any  case.  Inserting  the  value  of  b into  (3.12),  we 
get  the  minimum  value  of  X.  This  minimum  value  is  once  more  an 


<2  >.*rt  2 

increasing  function  of  (~~~  ) ** . The  reason  for  the  existence  of 

these  relations  is  as  follows:  In  a convective  fluid  motion  the 

temperature  gradient  is  a destabilizing  factor,  while  viscosity  and 
Coriolis  force  are  stabilizing  factors.  The  stabilizing  action  of 
the  Coriolis  force  is  studied  and  appreciated  in  dynamical  meteo- 
rology. In  (3.12),  the  second  term  of  the  right-hand  side  of  the 


equation  denotes  this  stabilizing  action  of  the  Coriolis  force# 

The  first  term  of  the  right-hand  side  denotes  also  the  stabilizing 
action  of  the  viscosity,  while  the  left-hand  side  of  (3.12)  shows 
the  destabilizing  action  of  inhomogeneous  heating.  Equation  (3.12) 
shows  that  the  convective  fluid  motion  can  exist  stationarily  when 
these  stabilizing  and  unstabilizing  actions  balance  each  other. 

OTaV-v^  o 

Thus  we  can  understand  the  \ (minimum)- — •( — relation  obtained 
above.  Furthermore,  as  is  shown  in  dynamical  meteorology,  the  stab- 
ilizing action  of  the  Coriolis  force  is  more  powerful  for  larger- 
scale  fluid  motions,  while  that  of  the  viscosity  is  more  powerful 
for  smaller-scele  motions.  These  circumstances  are  reflected  in 

(3.12) .  Thus,  the  first  (second)  term  of  the  right-hand  side  of 

(3.12)  becomes  larger  for  larger  (smaller)  values  of  b,  that  is, 
for  smaller  (larger)  3cale  motion.  In  short,  the  action  of  the 
Coriolis  force  makes  fluid  motions  of  a smaller  scale  easier  to 


excite  relative  to  larger-scale  motions.  This  is  the  reason  for 
- 2 

the  b ■ — ' relation  obtained  above.  From  the  above  physical 

considerations,  we  may  safely  assume  that  the  similar  relations,, 
namely,  X(min)  ^ end  b — '{ win  algp  be  obtained  ^adj^bci^dary 


conditions  other  than  those  assumed  here.  The  existence  of  the 

_ 2 

relation  b is  favorable  for  our  present  purpose.  The 


reason  Is  as  follows:  Since  the  Sg  -type  fluid  motion  Is  of  a 

smaller  scale  than  some  of  the  more  symmetrical  motions,  we  cannot 

expect  that  the  Sg^-type  motion  is  the  easiest  motion  to  be 

excited  (as  was  expected  in  ( 2 ) in  ^ 2)  without  the  relation 

b - — in  short,  while  the  S^0  motion  is  not  the  easiest 

motion  to  be  excited  in  the  case  when  (a  = 0,  we  may  expect  it 

2 

will  become  so  for  a certain  value  of  (~^*~)^<»  This  is  the  most 


important  result  obtained  in  this  section 


4.  We  shall  now  consider  convective  fluid  motions  in  a spherical 
vessel  of  radius  a.  Referring  to  the  rotating  coordinates  in  Pig.  2 
we  have  the  following  equations: 


A 


Fig.  2 

f>  ” grad  p + \iV2n  - Qp  (S  x u)  -a  g . 

M = k^2T  +@  ,(4.1) 

i ^ ^ 

where  u,  co  and  g denote  the  displacement  (of  the  fluid  relative 
to  the  above  coordinates),  rotation  (of  the  spherical  vessel)  and 
gravity  vector  respectively.  A source  of  heat  is  assumed.  The 
rate  of  he  it  generation  is  assumed  such  that,  in  the  absence  of 
heat  conduction  and  convection,  the  temperature  would  rise  at  a 
uniform  rate  0)  . Thus  in  the  state  of  no  convection,  we  have 

0 ~ = k^72T  + 

M ° 


In  a sphere  of  radius  a,  this  gives 


, ( 4 .0  ) 


where 


2\ 
r ) 


= JL 

2a 


( a 


(3  = 


HE 


(4.4) 


is  the  temperature  gradient  at  r = a.  In  (4.4),  Q and  {O ck  are 
the  heat  generation  (per  unit  time  and  volume)  and  heat  conduc- 
tiviyt,  respectively.  In  the  same  way  as  in  $ o,  we  get  from 
(4.1)  - (4.4) 


1 . 

v V u - -r  grad  p 


2(u»  x u)  + a ^T  = 0 


o dT 

v T = it  a?2  ‘ 


--k 


ru 
a r 


,(4.5) 

,(4.6) 


div  u = 0 


.(4.7) 


In  (4.6),  means  the  radial  velocity  in  the  spherical  co= 
ordinates  in  Fib.  2.  In  these  spherical  coordinates,  we  have 
also 

_ _ 

u = (u^,  UQ,  U^),  (a  = (a)  cos  <S,  - id  sin6>,  0), 

g = (g,  0,  0)  .(4.8) 

The  forms  of  velocity  vectors  u satisfying  (4.7)  are  known 
(H.  Takeuchi,  1150): 


u = 

(n  F 

+ r2d„ 

r 

' n,m 

XX 

r n&' 

\ = p 

f a_ 

j 0& 

r u4> 

J n,m 

9 

\ sin  e 

0) 


r ,m 


,(4.9) 
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where  a),  is  a solid  spherical  harmonic  of  degree  n and  is 

II 

expressed  as  follows 


co  = rns  (©,<£)=  rVH  6 Je1*^  .(4.10) 

n,m  n,m'  * r n 


In  (4„10),  s is  a surface  spherical  harmonic,  and  P®'  is  a 
n,m  n 


T?  — rs  4“  V*.  £1  TT/-\1  TT  TT  r\  A 4"  An  ~f*  4-  TT T*\  A 

* m i ' WU.  U j ■<.  | 


Legendre  function.  * or  the  velocity  vect 

dP dG 


1 wo  VioirA  ol  on 


div  u = 0 = 


n «.  n.m 


zS  + r . j*zH  + ( n+3 ) G. 


]» 

v*i  w J m t 


(4.11) 


(V2u). 


/l-\ 

/ a .x  \ 


/ . — ,2~>\ 

(V  u). 


a 

15z  / 


i 


/X\ 

(V  u)  - fn>m^r)|  9 ] “n,m  + gn,m  ( 7 ) Wn,m 

/ \ 2 / 


( v: 


r(V2u)  = (nf  + r^g  ) c»> 
xv  / u/  \ “-n.m  &n.m'  n. 


n,m 


3n,m 


.m 


r(V2u) 


e =.f 


/ 3 


'n,m 


‘30 

a 


< \sin  03?/ 


\ S ■ 
n»ai 

) X 


,(4.12) 


in  which 


n,m 


d_lFn,m,2n 
“dr"  ' r 


ni!£  4.  S> 


dr 


"n,m 


Sv 


2 

d ®n  ,m  j 2 ( n+2 ) d^n  ,m 

. a + — ? ax- 

dr 


.(4.13) 


Another  type  of  u satisfying  (4.7)  is  as  follows: 


Type  II 


n,ra 


/Ziy  “ yfe\ 

(r'  ' x3z  “ zax  | a’n,m 
.3  _3 


\y-5i 


3y 


-i: 


u 


tU  = 


*0 


U, 


0, 

3o 

T r‘ia_ 

ii,m  sin  db  <b 
Bo 

T n,m 

~ Ln,m  “5T” 


div  u = 0 


o 

t \ 7 £'„  > 

V V u x 

/yjsr  “ zl^\ 
— -*/  \ 

(V2u)  = 1 (r) 

v y n,m 

4-  - 4: 

3x  3z 

p 

s a / 

07  u)z 

yXB?  ~ / 

(V2u)„  = 0 

o 


n,m 


o Bo 

(V  u)»  = Wr)  sTn*tg3 


ou> 


/ r — 7 2 ^ n .m 

(V  u)  = I _ 53~ 

Y n.iu  ww 


d2L. 


n.m 


dr 


n.m  . 2(n+l)  ^n.m 
2 — ~v ar~ 


Taking  (4.5)  and  (4.6)  into  consideration,  we  may  assume 

T = T (r)o  , P = it  _(r)o  w 
n.mv  ' n.nr  n.m  n.m 

for  u of  type  I.  For  u of  type  II,  we  may  take 

T = 0,  F = 0 

The  Laplacian  of  T in  (4.18)  is  as  follows 


T7  T = t O , t 
v * n.m  n.ar  n.m 


a Tn.m  . 5(n+l)_  dln.] 

~~~2~  r dr 

dr 


Accordingly  equation  (4.6)  is  satisfied  by  putting 


(4.14 

(4.15) 


(4.16) 

(4.17) 

(4.18) 
, ( 4 • 19  ) 

1(4.20) 
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t 

n.m 


P_ 

ka 


(nP 


n.m 


r2G  ) 
n.m 


.(4.21) 


Furthermore,  if  we  omit  the  term  2(co  x u)  in 
fy  equation  (4.5)  by  putting 


v(nf  +■  r^g  ) 
' n.m  “n.m' 


Pv 


1 


n © tu 


(4.5),  we 


+ agrT 


n.m 


may  s at i s - 


= 0,(4. 22) 


Vf 


n.m 


1 

P ^n.m 


0 


,(4.23) 


1 = 0 or  L =0  .(4.24) 

n.m  n.m 

Thus,  when  co  = 0,  we  get  (4.11)  and  (4.21)  - (4.24)  from  (4.5)  - 
(4.7).  These  are  five  relations  among  five  unknowns  F , 

XI  elil 

G , L , T , and  it  • Vectors  of  type  I and  II  are  usually 
n.m*  n.nr  n.mJ  n.m 

called  poloidal  and  toroidal,  respectively.  Thus,  as  is  seen 
from  (4.24)  when  co"=  0,  we  can  satisfy  (4.5)  - (4.7)  by  a single 
poloidal  u and  the  corresponding  T and  P.  This  will  not  be  the 
case  when  S'  ^ 0.  The  existence  of  the  term  2(co  x u)  makes  the 
latter  problem  more  complicated.  A method  of  dealing  with  this 
complication  will  be  given  In  the  next  section. 

5.  We  denote  u in  (4.9)  and  (4.14)  by  u ( I)  and  u (II), 

n.  • m xi  • iii 

respectively,  and  shall  try  to  satisfy  (4.5). ? (4.7)  t)y  putting 

u = 21  u ( D + u ( H)  ..(5.1) 

*■ — n.m  n.m 

n 

*—  ^ 

The  d>  parts  of  ail  u (I)  and  u (II)  may  be  taken  to  be  cos  m 
T n.m  n.m 

(m  = 0,  1,  2,  . . .).  The  reason  why  we  may  assume  these  forms 
can  be  understood  from  (4.5)  - (4.20).  Equation  (4.11)  and  (4.21 
still  hold  for  each  pair  of  (Fn  m,  Gn  m,  Tn  m)  thus  introduced. 

We  now  make  the  following  transformations  of  variables 


-13- 


•(  = P (f  ) = V F l1')®  G if ) = G (r), 

' as  n.m  ’ 2 n,m'  ' s n.m  5 ' n.m'  5 

s. 

L (f  ) = L ( r- ) , T (f  ) = T ( r ) , u (f  ) = it  (a* > . ( 5.2  ) 

n.m'3  ' n.m'  ' * n.m1-3  ' n.m'  ' * n.m'  n.m  ' ' 


We  shall  denote 


cTF  (I)  ~ dF 

denote Eig—  - + -iySL  + 2G  -(f), by 

_,f  | a|  n.m 

5 

By  using  these  notations,  (4.11)  and  (4.21)  can  be 


by  f ( f ) 

^ n.m  J • - 


rewritten  as  follows 


t (f)  = - [" nF  (f)  +^2G  (/)?  (5.3) 

n.m'  3 k [_  n.m  3 J n.m  ' J 

d / c n v 

Fnf  ( f ) + £ 2g  ( / ) 1 — — -r  + £&£  g T (f)  = 0 (5.4) 

In  order  to  get  more  relations  among  unknown  functions,  operating 

with  i?  ( T)  (which  will  be  denoted  by  u^(  I)  hereafter)  upon  (4.5) 

'n.m'  p 

and  integrating,  we  get 

a n 2%  ? 

( ( ( (4.5)  • ur(I)  r sin  0 ard0d<£  = 0 .*(5*5) 

Jo  J o Jo  P 


'o  J o J o 


Inserting  (5.1)  and  (5.2)  into  (4.5)  in  (5.5)  and  carrying  out  the 


integration,  we  have 


5*  C(f  ) fZn?  (npPp  + t\)  at 


+ D(f  ) f 2ne  Fp  df 


,(5.6) 


n6 

C(i)  = (rj  f + f ) - i G(  1 1L(1 2 + 5&S  c T 

u^Xp+5  gp;  |hp-1  df  ’ v S 


n -n. 


- 2|§1  fl  (P)”"1  H | a p+1  xx  (P*a)La  =0  ,(5.7) 

L J a 


»<§>  = [VnP+1)  V^J 
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^ — n_-n„+l 

2_  I p = 0 


,(5.8) 


where 


I ( P ) = f (P«)a  sin  9 d© 


{ 5.9 ) 


i-^a-P)  ^ 


at 


s:  $ 


sin2©  a© 
dP  dP, 


(5.10) 


IgCa.p)  = I2(p.a)  * j paPp  + “TO  “S)sine  0030  d0  ^(5#I1-) 


m 


*“  (6),  f„  = F 


a n 


a 


rippiripj  — — 


(5.12) 


Similarly,  operating  with.  u^(  II ) upon  (4.5)  and  proceeding  in  the 
same  way  as  above,,  we  set 


,2(na+l) 

Lft  df-  0 


or 


F(?v  = 


{ *(f)f 

Jo 


2 

) = 1 -j.  2ga  r ^ J.T  'I  T ( Q\1 

1 -p  ■ v L“fj*--p 


(5.13) 


-1 


n -nft  *“1 

OC  p 


a 


[ Pa  + (naPa  + ^Ga}  I2(  a#p)j  .(5.14) 

If  equations  (5.3),  (5.4),  (5.7),  (5.8)  and  (5.14)  are  written 
down  for  all  P's  concerned,  we  see  that  we  have  sufficient  numbers 
of  equations  to  determine  the  unknown  radial  functions.  As  was 
said  immediately  after  (6.1),  the  values  of  m^  of  these  radial 
functions  are  constants  equal  to  a fixed  value  m,  say.  Furthermore, 
as  is  seen  from  (5.7),  (5.8)  and  (5.14),  the  velocity  vectors  u^ 
and  Uq  with  I_  ( a,  P ) , Io(a,p)  ^ 0 may  be  considered  to  be  coupled 
by  tho  action  of  Coriolis  force.  In  view  of  these  circumstances, 
we  may  call  I^(a.P)  and  3^(g,P)  the  coupling  integrals.  Now  it 
can  be  easily  shown  that 


15- 


? (nR+m R)  I 
j ( (3 ) - __iL — £_ 

ov  ' 2np+l  ( rip  "nip  )T 


(5.15) 


IgCa^p)  - + np(hp+l)  j PaPp  sinQcos©  d0  (5.16) 


and  that  the  coupling  Integrals  I^(asp)  ana  I?  ( a , p ) vanish  unles! 
nQ  - n^  = +1.  The  values  of  I^(a,p)  and  Ig(a,P)  for  several  a 
ana  p are  shown  In  Table  I.  Having  thus  established  the  funda- 
mental equations,  our  next  taSk  is  to  consider  the  boundary  con- 
'itions.  By  the  assumption  of  a "viscous”  fluid,  we  have 


u (I),  U (II)  = 0 at  f = 1 .(5.17) 

n.m  n.m  ' 5 

As  the  thermal  boundary  condition  we  shall  take  tentatively 

T = 0 at  | = 1 .(5.18) 

By  (5.18)  is  meant  that  the  temperature  at  r = a is  kept  constant. 
This  is  a plausible  condition  for  the  boundary  of  the  earth’s  core. 
Taking  (4.9),  (4.14)  and  (4.18)  into  consideration,  we  have 


nP  + % 2G  . 
n,m  n$ 


— u- 

m n .m 


L = T = 0 at 
n,m  n,m 


1 .(5.19) 


F , G 

n,m  n,nr 


- - must  also  be  finite  at  r = 0 (or  % = 0 ) . Thus 
our  problem  is  reduced  to  solving  (5.3),  (5.4).  (5.7),  (5.8)  and 


(5.14)  under  the  boundary  condtions  (5.19).  It  is  easily  seen 

a6g  a^ 

that  this  is  an  eigen-value  problem  for  A = o which  depends 

to  2 

on  the  parameter 


6.  By  using  the  results  of  the  last  section,  we  find  sets  of 

fluid  motions  as  shown  below 


U1.0* — 

— “2. O'1  ■L±^  ^ a3 .0^  ^ 

,,(6.1) 

~u2,0!  IJ  — ="*u5.0'  XI' 

,(6.2) 

U-1  t ( I ) < — 

X a JL. 

-*UP  t(II)^h=±u^  x(I)  - - - 

, (6.3) 
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Ulel(Il)  = 


■U2.1^  ^ 


-U3.1(  ^ 


,(6.4) 


U2  • 2 ^ X) 


-u„  Q(I±)«— — u„  0(  I) 

Omd  *±mf^ 


, ( 6 • 5 ) 


u2 .2 ' 11 


'U 


4.2 


(II) 


,(0.0/ 


u3.3^ ^ 


ru4.3(ll); 


5.3 


(I) 


,(6.7) 

etc. 


Generally  speaking,  there  exist  sets  of  fluid  motions  which  con- 
tain u ( I,  or  II)  with  n = m as  their  first  member.  As  is 

n •in.  * 

easily  understood,  the  case  when  n = m = 0 is  an  exceptional  one. 

This  is  the  most  important  result  obtained  in  the  last  section. 

2 2 

0(  I)  in  (3.5)  is  usually  denoted  by  S0  and.  it  is  this  S0- 

type  velocity  which  was  referred  to  in  section  2.  In  a rotating 
2 

sphere,  the  S0-type  motion  cannot  exist  by  itself.  It  must  be 

accompanied  by  the  fluid  motions  of  1S3  g(  II)  , g(  I),  - - - 

2 

types.  We  shall  call  these  latter  quasi-Sg-type  motions.  For 
2 

the  quasi-Sg-type  motion,  we  have  from  (5.3)  and  (5.4) 

dG, 


2 d^2  '"'"2 

j ar  + f ar  + “2  = 0 


4 ®4 


dG, 


f Tf  + f Hf  * W*  = 0 


d T 


df 


d2T 


a dT  p p - d r, 

T-  * f Tf  = - <2p2+t“G2) 


dT 

jr  4.  " ■* 

TW 


4 A ID  ai4  _ Paw  ^ 

rxc — — (^ 


/°4> 


.4(4) 

,(b)  etc. 

(6.8) 

, ( a) 

,(b)  etc. 
( 6.9) 


In  (6.8)  and  (6,9),  we  write  simply  FP, 
Jp  o< 


or 


2.2 


(#) 


ii'  2 

Go  o(f),  “ - -.  Corresponding  to  C(f ) = 0 in  (5,7)  we  have 


r d2Pp  . dF  2 d2G, 

! 2{Tf-  + ?TT  + 2a2}  + ^ 


w 


df 


2 , 8 dG2  _ 1 d^  ^2  ^ 
I ~ df  * 
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. agoa  e 2 5 n , 192.  £-2.  _ 

+ / T2  + 23  J L 1 T->  ? L3  - 


0 


(a) 


4 ( _ •*■  + 2.  /*  + 2^  ' + f —i  + i—  __i)  1 4 

4(C5F  s *r  4'  § 


d2p4  „ 8 *4 


pvf 


~dT 


a®oa  - 2_  1 n /4 80t  480x7  c2t  ^ ^ / v,  x „ „ 

+ ~v — 5 T4  ” IS  (_T^3 H I L5-'  ~ ° ,(b)  etC# 

(6.10) 


In 


(6.10)  SI  = ^ 


irid  g is  put  equal  to 


S = g„f 


(6.11) 


,(6.12) 


where  gQ  is  the  gravity  at  f — 1,  that  is,  at  the  surface  of  the 
earth's  cor9.  The  gravity  g in  a uniform  self -gravitating  sphere 
may  be  shown  to  vary  as  in  (6.12)-  Next,  corresponding  to 
D(f)  = 0 in  (5.8),  we  have 


2 

d F, 


dS‘ 


S+i.^k+20  - !£  - -i.il  22if2L  =n  (a) 


dSp4  + 8 *4  . 20 

+ _ + 2G4  - — STO 


*4  - ( .4.0Q?-5-L5  + i§g28  |2L5)  - o ,(b)  - etc. 

(6.13) 


Lastly,  corresponding  to  E(  J ) = 0 in  (5.14)  we  have 

L3  . 8 dL3  . 7 _n  (%-2ir  192  . 384  . 2p  +fB  J 

ST  + J H|“  + I44U  J L i “7“  F2  + ”T“  UF2  + 5 VJ 

rf i«4  - $ u -4ij 


d2L. 


480  F4  + (4F4  +1^), 


= 0 


,(a) 


i2jl5  12  dL5  . 11 

"Ts-  + y 3T 


af 


.12:1  (s-y.  ^ + 480x28 


. T 2800x24  „ . 2800x168  ,m  A „2„  , 

+ T Hir  p6  + ”“11x13"  voP6  + 1 V 


/ = 0 - (b ) etc. 
) (6.14) 


If  we  omit  u 0(  I,  or  II)  with  n ^ 4 in  (6.5),  we  get  the  results 
(a)  in  (6,8)  - (6.14).  In  that  case,  the  F^  and  G^  terms  In  (6.14 
(a)  must  of  course  be  omitted.  The  results  thu3  obtained  will  be 
called  the  approximation  (a).  Approximations  (b)  and  (a)  will  be 
defined  similarly*  In  the  approximation  ( a)  we  have  five  unknowns 
Fg,  G^,  -jTg,  Tgjj  Lj.  In  the  approximation  (b)  the  number  of 
unknowns  becomes  10, 

proceeding  in  this  way  to  approximations  (c)  and  ( d ) -v.  * - w 
shall  obtain  the  required  solution  for  (6.5).  Boundary  condi- 
tions to  be  satisfied  are  shewn  in  (5.13).  Although  the  above 

2 

results  are  obtained  for  the  quasi-S0-type  motion,  similar 
results  may  be  obtained  for  (6.1)  - (6.7).  It  is,  however,  almosi 
hopeless  to  solve  the  eigen-value  problem  by  trial-and-error 
methods.  A method  of  dealing  with  the  difficulties  will  be 
given  In  the  next  section. 

7.  Taking  (5.19)  into  consideration,  we  shall  put 


where  A , 
o' 

ing  (5.5) 


"hum  +^an.m  = (1-5>2  <*o  + *1#  + >> 

Lium  = ’(W)  (B„  +Bif+  - - -) 

B , A. » Bn  - - - are  undetermined  constants, 
o 1 1 

into 


( n+1  }G 


d(nF 


n.xn 


+ f G ) 
n.m 


n.m 


TTf 


•(7.1) 

Transform 


,(7.2) 


we  can  calculate  the  value  of  G by  (7.1)  and  (7.2).  Inserting 

n •m 

G thus  obtained  Into  (7.1),  we  get  Y . The  form  of  nP  + 
n.m  \ o n.m  n.m- 

O 

in  (7.1)  is  so  adjusted  as  to  make  F vanish  at  § - 1. 

" u.xii  n.m 

Next,  inserting  (7.1)  into  (5.4),  we  have  the  following  differ- 
ential equation  for  T 

n.m 


-19- 


d2T  o,  .-1X  d.T  ..3 

n.m  , 2(ntl)  n.m  po, 


df‘ 


df 


1C 


V Ar 


,(7 .3) 


'.i 


sphere  A*  is  a certain  constant  which  is  a function  of  AQ,  - - - 
in  (7.1).  The  equation  (7.3)  is  solved  in  the  form 


Tn.m(1ir)’i  “ A"  " ■(j+S-)Td+3+5HT 


,(7.4) 


...1 a tl  A ~ 

WiitJX'c?  a a.  3 


a constant  which  is  determined  so  as  to  make  T in 

n.m 


(7.4)  satisfy  the  condition  (5.19).  The  value  of  A"  thus  adjusted  is 


Tj+S')(  j+3+gn)' 


Inserting  (7.1)  and  (7.2)  into  (5.8),  we  can  calculate  it 


n.m 


(7.5) 

Thus, 


with  nF  4=  ^ and  L in  (7.1).  to©  can  satisfy  three  grouns 

n.m  n.m  n.m 

of  differential  equations,  (5.3),  (5=4)  and  (5.8).  There  remain, 
however,  two  groups  of  differential  equations  to  be  solved.  They 
are  (5.7)  and  (5.14).  We  have  also  two  groups  of  undetermined 
constants  AQ,  A^$  ~ “ -Bq,  - - -.  We  shall  now  try  to  satisfy 
(5.7)  and  (5.14)  approximately  by  choosing  these  constants  ade- 
quately. Inserting  (7.1)  - (7.5)  into  (5.13)  and 

,(7.6) 


p ji , 0 

{3,5.  -r.  « . ^ 2 , 


j C(f  )§  P(fipFp  + f Gp)d|  = o 


(see  (5.6)  and  the  equation  D(£ ) = 0 above  satisfied)  we  have 
equations  of  tne  following  forms 


o \ 

'O 


C1  C(f  )f2nP(i-|)2df  = 0, 


2n„+l 


- -u  » a-  q 

A1  ) C(^)f  ^ U-§)~df  = °, 

A c 2(n  tl) 

Bo  f E(f)§  P (l-r)df  = 

- o 


,(7.7) 


0 
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rl  2(nfl+l)+l 

B-,  \ l£(§)f  1 (i-f)dj  = 0,  - - - .(7,8) 

Jo 

These  equations  will  be  used  to  determine  the  values  of  Aq,  A-^,  --- 

Bq,  B^,  - - -•  In  order  to  show  the  way  of  determining  AQ , - - 

we  shall  take  up  the  approximation  (a)  fbr  the  case  (6.5).  We 

2 

shall,  furthermore,  take  only  the  first  terms  of  2F0  0 + f ^2  2 
Lc  c in  (7.1).  The  equations  that  determine  A and  B ere  obtained 
by  (7.7)  and  (7.8),  as  follows; 


/ 571  -v  2 w 

•11,  18.  30.  Si.  SST  13  A “ 15;ac 


, 2 Pi 


-fl 


I, 


(ST  T5TTF^Ao  * (”  9)b< 


+ < §T7S)Bc 


o 


o, 


.(7.9) 


In  order  that  these  equations  be  compatible,  the  determinant 
formed  by  the  coefficients  of  A and  B must  be  equal  to  zero. 

A _ ^ w 

a^a(3g 

The  equation  for  'A  = . — ^ ..  thus  obtained  is 

^ _ 2 ’ ; ' • 

A = 10602.5  + 1 .86994 12  ^ , iZ  = 1(7.10) 

It  is  to  be  noted  that  this  is  a relation  of  the  same  type  as 
( 5 • 12  ) • The  values  of  A in  (7.10)  have  been  calculated  for  several 

/'S  O 

j 's  and  are  shown  in  the  column  (a.l)  in  Table  II.  The  results 

V q 

obtained  by  taking  the  first  two  terms  of  2Fg  ~ + | Gg  g I»2  0 

are  shown  in  the  column  (a,2)  in  the  same  table.  Similarly,  the 
(b.l)  and  (c.l)  values  of  A are  shown  in  Table  II,  From  the 
results  in  this  table  we  see  that  we  m ay  u s g t ho  (c.l)  values 
of  A for  the  exact  A in  the  case  (6.5)..  The  (c.l)  values  of  A 
for  the  cases  (6.1)  - (6.7)  are  calculated  and  are  shown  in 
Table  III.  A discussion  of  the  results  given  in  Table  III  will 
appear  in  the  next  section. 

8.  As  is  seen  from  Table  III,  if  we  put  o’  = 0 in  (6.1),  (6.2), 

4 4 

- - - . we  get  the  eigen-values  A = 0,811  x 10  , 1.06  x 10  , - - 
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The  eigen-functions  corresponding  to  these  eigen-values  are  of 
« o 1 1 2'  2 ^ 

pure  Sf,  3.-lS  Sn  , S~,  S0 , S~,  and  3-r  types.  These  are  nothing  but 

X c*  j_  {1  O O 

the  results  obtained  by  H.  Jeffreys,  M.  E.  M.  Bland  and  S.  Chandra- 
sekhar ( H.  Jeffreys  and  M.  E.  M.  Bland,  19  51:  S.  Chandrasekhar, 

s~\  2 2 3 

1952).  The  minimum  values  of  A in  the  cases  J L = 0,  10  , 10v  and 

104  are  given  by  the  fluid  motions  of  the  quasi-8^,  S^,  S^,  and  Sg 

types,  respectively.  The  minimum  A increases  with/2.2.  These 

2 


relations  correspond  to  the  b 


2csh2v2 


)CJ  and  A(min)-~^  (“-— ) relations 


of  ^ 3.  From  our  present  point  of  view,  it  is  important  that  the 
2 

quasi-Sg-type  motion  becomes  the  easiest  type  to  be  excited  when 


n 


L * = 10 


4 


We  shall  now  estimate  the  value  of  the  kinematical 


viscosity  v corresponding  to 


a 2 = io4,  n f *£  - 1°2 


(0.1) 


Putting  (0  = Q-g^QQ  (sec)”'*'  and  a = 3.4  x 10°  cm  into  (8.1),  we  have 
"112 

v = ID  cm  /sec.  In  short,  if  the  kinematical  viscosity  of  the 

11  2 

fluid  in  the  earth's  core  is  about  10  " cm  /sec,  we  shall  have  the 

2 

fluid  motion  of  the  quasi-Sg-type . This  value  of  v Is  very  large 
compared  with  that  which  was  hitherto  estimated  for  the  earth's  core, 

i.e.,  v = 3 0 1 — '10  ^ cm2/sec  (E.  C.  Bullard,  1949).  This  contradic- 

-2  -3 

tlon  may  be  avoided  as  follows.  While  the  viscosity  v = 10  ^—'10 

above  referred  to  is  the  molecular  one,  what  our  present  study  is 
concerned  with  may  be  some  kind  of  equivalent  viscosity.  In  our 
present  study,  for  example,  the  existence  of  a magnetic  field  in 
the  earth's  core  and  its  prohibitive  action  on  the  convection 
(S.  Chandrasekhar,  1952  ) have  not  been  taken  Into  account.  The 
equivalent  viscosity  of  the  magnetic  field  may  be  estimated  as 
follows:  The  equations  of  motion  in  a magnetic  field  are 


an  _ 1 

“ Wfi 


( curl  K x H)  + 


vV2u 


+ - - - 


,(8.2) 
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The  assumption 


4^3  ( curl  H x H)-~  vV^u 


(8.3  ) 


may  be  used  to  obtain  an  expression  for  the  equivalent  viscosity 
ox  the  magnetic  field, 

- ^ .(8,4) 


V ~ 45^5 

In  (8.4),  L,  H ard  U are  representative  values  of  length,  magnet it 
field  and  velocity.  Inserting  - 10,  L = 3 x 10^,  H = 40  and 
TT  = 0.03  into  (8.4),  we  get  v^.bxS.C^'  cm^/sec  which  is  of  the  orde 
of  magnitude  required  in  (8.1).  We  shall  now  consider  the 
corresponding  value  of 


\ = 2.3  x 104 

obtained  in  the  last  section.  Inserting  gQ 
a - 10*"5,.k  = 10"1  and  v = 1014  into  ^8.5),  we  get 


(8.5) 

103,  a = 3.4  x 108, 


p = 2 x 10”16  -£■  (8.6) 

r cm 

As  was  stated  in  section  3,  p is  the  temper ature  gradient 

(at  the  surface  of  the  earth's  core)  in  excess  of  the  adiabatic 

_5 

one.  The  latter  is  estimated  In  the  core  to  be  10  Thus  the 

O 

quasi-Sg-motion  can  exist  In  a stationary  state  with  a thermal 
gradient  which,  differs  but  little  from  the  adiabatic  one.  In 
estimating  p in  (8.6),  we  used  the  value  of  the  molecular  thermal 
diffusivity  for  k.  Even  if  we  substitute  some  kind  of  eddy 
diffuaivity  for  it,  the  res"it  obtained  above  will  not  be  changed 


much 
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Table  11  A in  lu  . 

PI* 

(a.l) 

(a. 2) 

(b.l) 

( c.l) 

0 

1.1.060 

1.06 

102 

1,079 

1.08 

1.078 

1.078 

1 o5 

1.247 

1.27 

1.236 

1.235 

o 

H 

2.930 

3.03 

2.359 

2,311 

Table  III 

X in  104 

Pi 

ui,o( 

U1.0^ II' 

ul.lU) 

u^idl) 

U2.2^  -1' 

U2.2' 11 

U3.3(  1 
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1.06 
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1.06 
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1.53 
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0.862 
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0.849 
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3.52 
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